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Using effective field theory techniques we calculate the source multipole moments needed 
^\ ■ to obtain the spin contributions to the power radiated in gravitational waves from inspiralling 

compact binaries to third Post- Newtonian order (3PN). The multipoles depend linearly and 
quadratically on the spins and include both spin(l)spin(2) and spin(l)spin(l) components. 
The results in this paper provide the last missing ingredient required to determine the phase 
evolution to 3PN including all spin effects which we will report in a separate paper. 
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INTRODUCTION 



Coalescing compact binary systems are one of the most promising classes of sources in the search 
for gravitational waves through direct detection at earth based interferometers such as LIGO [l|], 
VIRGO [2j and at the planned LISA [3j and Einstein Telescope 4] experiments. One of the main 
goals of these instruments, beyond wave detection itself, is to extract masses and spins from the 
signals. However, the ability to perform these extractions is predicated upon precise theoretical 
predictions. 

During the initial inspiral phase of the binary coalescence the compact objects are well separated 
and one can study their dynamics analytically in the post-Newtonian (PN) approximation. For 
binary systems the PN approximation is a non-relativistic expansion in the typical three- velocities 
v in the binary, and the standard lore is that accurate parameter extraction requires the knowledge 
of the phase of the gravitational wave signal at least to 3PN 1 or 0(v®) [5rtZ]. For spinless systems 
3.5PN precision has been reached in the last decade by a series of dedicated calculations, see 
for a review and further references. 



Recent measurements of black hole spins indicate that black holes in binary systems may fre- 



quently be close to maximally rotating [9], and that spin effects have a significant impact on the 



development of accurate gravitational wave templates 



10M13I] . In addition, higher order PN correc- 



tions including spin are relevant for comparison between analytic results and numerical simulations 
[14J . Il5j ] . It is therefore desirable to account for spin effects at the same degree of accuracy as the 
spin-independent components. 

Currently, the phase evolution of the expected gravitational wave signal is known to 2.5PN 



precision when including spin (in a power counting assuming maximally rotating black ho 



This includes the leading order (LO) spin-orbit terms, i.e. effects 
the LO spin(l)spin(2) and spin(l)spin(l) contributions at 2PN 



inear in spin, at 1.5PN 



order (NLO) spin-orbit effects at 2.5PN 



m]2M- 



16 



es 



12|, 



16 18] and the next-to- leading 



1 We follow the usual convention in the literature and label leading effects as OPN with higher order corrections 
counted relative to these. For example, the leading conservative dynamics gives the Newtonian or OPN equations 
of motion, the leading energy flux in gravitational waves due to the leading spinless mass quadrupole (OPN) of 
the system is called the OPN flux, and the phase evolution which follows from the OPN conservative dynamics and 
flux is referred to as the OPN phase. 

2 The spin S of a compact object is related to its mass m by |S| = aGm 2 where G is Newton's constant and < a < 1. 
Assuming maximal rotation means we power count a ~ f , but our results hold of course for arbitrary a. For binary 
systems in the PN regime and for our purposes the spin scaling is conveniently written as |S| ~ mrv 2 ~ Lv where 
r is the orbital separation between the binary constituents and L the orbital angular momentum of the binary. 



In this paper we calculate all missing ingredients to obtain the spin contributions to the energy 
flux and the phase evolution to 3PN order using the effective field theory (EFT) techniques first 



introduced in 



22, 



231 ] . and extended in |24|, 125J] to include spin degrees of freedom. We include 
the following multipole moments: The mass quadrupole with 0(Sa) and 0(S A ) components to 
NLO as well as the leading 0(SaSb) terms. The current quadrupole with LO 0(S A ) and up to 
NLO 0(Sa) contributions. The 0(Sa) and 0(S A ) terms for the mass octupole at LO. The LO 
effects of 0(Sa) in the current octupole, and the spin-independent 1PN corrections to the current 
quadrupole moment. 

While the c omp onents of 0(Sa) have been computed at the same order but with different 



conventions in 



2ll j. the NLO 0(S A ) as well as the 0(S^Sb) terms are new results. Together 



2614301] using the same conventions as in 



with the conservative part of the dynamics obtained in 

this work, our results enable us to compute all spin-dependent components of the energy flux and 

thus the phase up to 3PN. These will be reported in a separate publication 



include a comparison with the 2.5PN results of 



2l| 



311 ] where we will also 



Our computation is performed within the framework of Non-Relativistic General Relativity 
(NRGR) [22]. The EFT treatment (see 32] for a review of EFT methods) of binary inspirals is 
based on the separation of the physical scales in the binary which are disentangled and treated one 
at a time. This is achieved by introducing a nontrivial decomposition of the gravitational field into 
different kinematic regions and by setting up several EFTs in stages removing each scale. This 
makes NRGR a systematic and simple tool to perform calculations in the PN expansion, see 
for a pedagogical review. 



33] 



The original paper proposing NRGR 22| contained the computation of the spinless 1PN conser- 
vative dynamics and the LO radiated power which is due to LO spinless mass quadrupole radiation. 



Dissipative effects were included in 



341 ] via the introduction of additional worldline degrees of free- 



dom. The EFT framework works very efficiently in higher order PN calculations and it was used 
351 ] to compute the spinless conservative part of the binary dynamics to 2PN or NNLO using 



m 



the convenient parameterization of the gravitational field of J36l ] 



371 ] where radiative corrections 



The EFT for the radiation sector has been developed further in 
to gravitational wave emission such as tail effects were considered including a proper treatment of 
all divergences. The power emitted in gravitational waves from spinless binaries was reproduced 
at 1PN, 1.5PN and 2.5PN, and the leading ultraviolet logarithms in the gravitational wave power 



which arise at 3PN were computed and resummed using the renormalization group [37l ] . In order to 
calculate real-time quantities such as the gravitational waveform or the self-force in a self-consistent 
framework, NRGR was translated to the in-in formalism in [38 . 

NRGR was extended to include spin degrees of freedom in [24[, where the LO spin-spin and spin- 
orbit effects were reproduced as a sum of very simple Feynman diagrams. The NLO spin(l)spin(2) 
potential was then first derived u sing the Newton- Wigner (NW) spin supplementarity condition 



(SSC) at the level of the action in 



2514271] . Afterwards in [28] these results were re-derived using a 



Routhian approach in the covariant SSC, which was also used to compute the NLO spin (1) spin (T) 



potential including finite size effects 



29]. These results were independently corroborated in 39, 4C] 



using a more traditional set of technicmes, and complete agreement on the conservative 3PN spin- 



spin dynamics has been reported in 4l( . The NLO spin-orbit potential was recently derived using 



the EFT approach in 



30, 



42 



431 ] and shown to agree [30| with the results in [20J . Dissipative effects 



due to spin were also computed in 



i|. 



The worldline EFT techniques of 22| are clearly not limited to the PN regime of binary 
inspirals, or only to gravitational interactions, and can be applied to other systems with separate 
scales where a sensible power counting can be established. Inspired by 12211 . these techniques have 



been used to study the electromagnetic self-force on extended objects 45J], the radiation-reaction 

46 ] and the thermodynamics of caged black holes [47H49]. 
More recently, an EFT setup has been proposed to study cosmological perturbations [501 ]. 



force in the extreme mass ratio limit 



This paper is organized as follows. In Sec. [IT] we summarize the EFT setup for radiation effects 
from spinning extended objects. The matching calculation for the source multipole moments is 
performed in Sec. IIIII Our combined results for the multipole moments are presented in Sec. IIVI 
where we outline the assembly of all ingredients to the 3PN spin effects in the radiated power. 
The main results of this paper are collected in Eqs. ([791 - [86]) written in the covariant SSC. The 
energy flux and phase evolution will be reported in a separate publication [311 ] . We conclude in Sec 
IV1 The Appendices contain some subtleties regarding spin effects in GR and the spin-independent 
multipole moments computed with our set of conventions needed for the power loss to 3PN. 



II. EFFECTIVE FIELD THEORY SETUP 



The physics of a compact binary system is governed by three relevant scales, the size of the 
compact objects R, their separation r and the wavelength A of the gravitational waves it emits. In 



the PN regime all three are well separated and their hierarchy R <C r <C A is controlled by powers 
of the small expansion parameter v. We have the relations R ~ rv 2 and A ~ r/v where the former 
assumes compact objects (such as black holes). 

The NRGR construction of EFTs describing binaries proceeds in stages. Above the scale R 
we describe the compact binary system as two point particles coupled to gravity with higher 
dimensional operators 3 encoding finite size effects. This EFT is valid in a general frame, and can 
be written down in a generally covariant fashion. The corresponding action can be utilized in either 
the extreme mass ratio or PN scenarios. 

In order to properly disentangle the scales we decompose the gravitational field into different 
momentum regions which are responsible for different physical effects. Hence h^ v jra v = g^ u — r]^ 
(where the Planck mass m p is related to Newton's constant by G = l/(327rrrip) in our conventions) 
is written as the sum of potential modes H^ u and radiation modes h^ u with momentum scalings 
daHfn, ~ (v/r,l/r)Hfj, u and dah^ ~ {v/r,v/r)h^ u respectively [22]]. The potential modes are 
responsible for the binding dynamics of the binary whereas the on-shell radiation modes describe 
the gravitational waves which propagate out to the detector. Moreover, the radiation field hn V has 
to be Taylor expanded around a single point in order to achieve a uniform power counting at the 
level of the action [5l|] . By working in the background field gauge one can preserve diffeomorphism 
invariance in the radiation sector of the theory, and furthermore the potential and radiation modes 
can be gauge fixed independently 4 . 

To compute radiation observables, such as the energy flux, we "integrate out" (i.e. solve for 
and plug their solution back into the action) the potential modes. The resulting radiation theory 
describes the binary as a single point particle coupled to gravity and endowed with a set multipole 
moments which describe its internal dynamics. 

In the next section we review the basic aspects of the EFT approach. Readers interested in 
further details should consult 



22 



24 



37] 



3 The term "operator" is used here only to be consistent with terminology developed for general EFTs. 

4 There are no double counting issues between modes in the EFT due to the differing analytic structure of their 
respective propagators in momentum space (see [32|] for details.) 




FIG. 1: Leading order contribution to the effective action after integrating out the radiation held. The 
double solid line denotes the binary system and the wiggly line denotes the radiation graviton. 

A. The Radiation Sector of NRGR 

The EFT describing the radiation couplings of the binary system is naturally formulated as a 
multipole expansion, since it corresponds to a long wavelength approximation. It is convenient to 
work in the center-of-mass frame where the action of the radiation EFT [37| reads 

S Ta d[h] =- dt( m\fm + -xL lj uj ij0 J 

+ \j dt f&Eij - p ij B i:j + l -I l i k V k E l3 - ^ k V k B tj + • • • V (1) 

The first two terms contain the couplings to the total mass m and the total angular momentum L tJ 
which do not radiate since they are conserved. However, these terms source the Kerr background 
due to the binary itself in which the emitted gravitational waves propagate. Note also that our 
definitions of mass m and angular momentum L lJ are equivalent to the ADM prescription. The 
source multipole moments which can radiate are of either electric or magnetic type parity, and 
couple to the electric (Eij) or magnetic (Bij) components of the Weyl tensor and derivatives 
thereof. 

The observable we are ultimately interested in, the power loss or energy flux in gravitational 
waves, can be either extracted from the imaginary part of 'vacuum energy' diagrams such as the one 
shown in Fig. [1] or from squaring graviton emission amplitudes and summing over polarizations. 
Either way yields the standard result 





5 \ V^ 3 J / 45 \\dt 3 J / 189 \\dt 4 v ' ) I 84 \\dt 4 

'(2) 

where we did not include effects in the radiation theory due to the nonlinearities of gravity such 

as the tail effect. 

The information about the binary and its internal dynamics is encoded in the Wilson coefficients 

of the action in Eq. ([!]), i.e. in the multipole moments. These are obtained through a matching 



calculation which is the central part of this work. For this purpose we compare the action for the 
radiation EFT in Eq. ([!]) with the EFT valid below the orbital scale r, but above the finite size scale 
R, which describes the binary as two point particles coupled to both potential and radiation modes 
of the gravitational field. We then calculate the one-graviton emission amplitude from Feynman 
diagrams with a single external radiation field and arbitrary potential exchanges between the two 
point particles. This amplitude can be expressed in the form of an effective action linear in the 
radiation field h^ v 



T[h] 



1 



2m r 



d A xT^{x)h^(x), 



(3) 



where d tl T^ v = due to the gravitational Ward identity. Eq. ([3]) serves as our definition of the 
stress energy pseudotensor T^ v which contains gravitational contributions due to potential modes. 
In order for this effective action to obey a manifest power counting the radiation field hp. v {x) 
must be Taylor expanded around a single point within the binary [51] which is conveniently chosen 
to be the center-of-mass x cm = 0. Plugging this Taylor expansion into Eq. ([3]) and organizing the 
multipole moments into irreducible representations of the rotation group using symmetric trace- 
free (STF) tensors (plus relations which follow from dJT^ = as well as the on-shell equations 
of motion), we can bring Eq. ([3]) into the same multipole expansion form as the action for the 
radiation EFT of Eq. (fTJ. This procedure yields the matching for the multipole moments needed in 
terms of moments of T^ v . At the order we are working, the following expressions for the multipole 
moments are sufficient 



1 o 



J? 



1 2 



i 



Ajk 



J 



r tJ 



J, 



tjk 



d 3 xT V* J ]TF, 



hOll 



11. 



00^2 



cFx T" - -T U V + — T uu x z x'x 
3 42 



A v r 



TF; 



ci 3 x f-T"x 2 + lf lm ^ m - -T ,C Vx 2 + — r w x 4> l [xV 
V21 6 7 1512 J L 



d 3 xT 00r x * x J x fc 



TF, 



d 3 x f T u - f ol K l + -f 00 x 2 



x i x J x fc 



TF 



d 3 x 



d 3 x 

d 3 x 



e i/fcjiOfc x J x « 



STF 



1 3 

ilk I x j'^x^'x^x" 1 4- — T^x'x 2 

28 



28 



STF 



TF 



STF 



(4) 
(5) 
(6) 
(7) 
(8) 
0) 
(10) 

(11) 



where the STF prescription only acts on the free indices of the multipole moments. Note that the 
couplings of the multipole moments which can radiate, quadrupole and higher, are all contained in 
the hij part of the effective action in Eq. (|3|) since these describe the physical on-shell modes, i.e. 
the gravitational waves. Thus, in obtaining Eqs. (fi llip we may Taylor expand only the radiation 
field hij and set Jiqq and hoi to zero. 

In computing the stress energy (pseudo)tensor which enters in Eqs. (jH- llip it is convenient to 
consider the partial Fourier transform 



T^(i,k)= / d 3 xT^(t,x)e 



-ikx 



(12) 



since the calculation of Feynman diagrams is most easily performed in momentum space. In 
particular, we can read oST^ v (t, k) by calculating all diagrams with one external radiation graviton 
leg. By formally Taylor expanding 



T^(t, k) = f; ^f- ( f d 3 xT^(i, x)x* • • • x^ k n • • • k ln 

71=0 n ' ^ ' 



(13) 



we can extract the moments needed in Eqs. (jH - [TTj) by direct comparison with the results from 
the diagram. For further details on the radiation EFT and the matching see [371 ] . 



B. Spin Couplings 

The main advantage of our method is the introduction of an action from the onset. A spinless 
point particle of mass m is described by the familiar worldline action 



£p P = — m I dr 



(14) 



with dr = y / g^ lu dxi J -dx u , and spin is added as an additional worldline degree of freedom. 



couplings were introduced in the form of an action in 



24, 26L Mora Routhian in 



M 



28, 



Spin 



29]. In 



both cases the leading spin-graviton coupling is written as 



Ss g = ~\j S ab oj abfl u»dt, 



(15) 



in terms of the spin in the locally flat frame S = S^ v e c t L e u and the spin connection uj a bn = e-vbD^e 1 ^, 



where e* is a vierbein such that e*e"r/ a fe = g^. We chose A 



t parameterizing the worldline 



10 



in terms of coordinate time, and therefore u a = ^- = v a = ^|- = (1, v). 

Now we expand Eq. (|15|) in the weak gravity limit. We note that after expanding in powers of 
hnu we do not distinguish between local and covariant indices any more, but we consistently work 



with spin in the locally flat frame. We obtain at linear order in h 



{"' 



L 



LI 



LI 



2m, 



hio,kS , 



oft 



- — (hi jjk S ik v j + h 00 ,kS' 
2m p \ 

2^- (h oj>k S° k vJ + h io ,oS io ) , 



LI 



sg 2m, 



■hijflSPv*, 



(16) 

(17) 
(18) 
(19) 



where we kept the counting of v for potential modes of 

S^yi) both scale as relative factors of v. At 0(h 2 ) we have the terms 



24) ] . Time derivatives or factors of S l0 (^ 



L 



Sh 2 



—S a0 h^ I -K\ jfl + h^a - V*,a J u^ 



4rrip 



(20) 



Splitting hnu hito the sum of radiation (hn V ) and potential {Hnu) niodes we obtain the Feynman 
rules linear in spin. The terms linear in the metric are depicted in Figs. [2jz for radiation and in[2f> 
for potential modes respectively. 




i) 



b) 



FIG. 2: 0(Sa) worldline coupling to: a) radiation and b) potential modes. A blob represents couplings 
linear in spin induced by the terms in L sg . 



1. Spin Supplementarity Condition (SSC) 



The choice of SSC plays an important role in our analysis. Since we work with spin in the local 
frame, the SSC inevitably involves the particle's velocity in the local frame, and this introduces 



important effects in the dynamics 



27 



281 ] since we wish to express the velocity in the PN frame. 
For instance, to the order we are working in this paper the covariant SSC S Pb = can be written 



11 



as S e&^u^ = 0, and we can express 



Da 



H 



ei(xA)v^ + e^(xA) 



1 + 2(1 



4M 



SJ + - 



(21) 



with k = 1 and 1/2 for the covariant SSC and Newton- Wigner (NW) SSC respectively. To the 
order we work in this paper the vierbeins in Eq. (|21[) are given by 5 



eg(xA) 



4(*a) 



GrriB 

r 

^ + 5(nxS s f + 



G 



4 M _(l + «2-)*, 



+ 



(22) 
(23) 

(24) 



with r = Xj4 — X£, r = |r| and n = r/r. 

Throughout the whole paper we impose the condition B ^ A whenever we employ the index B 
as a label for a binary constituent. With these we obtain 



rf? 



Kb* 



^ G (^„uJ,J^Lt^*^\^<^^±^ 



v J A + - 2m B (v J A - v J B ) + -(r x S B )^ + (1 - k 



+ 



(vxxS.OMl + fi-^U^ 



G 



((v>4 - v B ) x S A ) 1 + ^ (S J B r • Sa - r J S A • S B 



+ 



(25) 



for the on-shell relation between S^ and the spin 3-vector in the NW and covariant SSCs 
The alert reader may wonder at this point about the extra term in the Routhian of 25, 



28, 



m, 



schematically of the form RSSuu, needed to ensure the SSC is preserved in time. Even though this 
term vanishes on-shell, it does contribute to the equations of motion and to radiation. However, one 
can show that this extra piece does not enter in the multipole moments until 4PN, see Appendix 

El 

C. Finite Size Effects 

The LO spin(l)spin(l) contributions to the multipole moments arise due to finite size effects. 
Contrary to the case of conservative effects, the 0(S\) contributions at NLO needed for the 3PN 



The effects on the dynamics (and on the preservation of the SSC) due to radiation do not enter until higher orders. 



12 



energy flux arise purely due to finite size operators. Here, we give an overview of the necessary 



Feynman rules for these terms, see 29J for further details. 

In the EFT framework of NRGR finite size effects are naturally encoded by a set of higher 
dimensional operators in the worldline action which can be written in terms of the electric and 



22 



24 



341 ] . In addition, these operators 



magnetic components of the Weyl (or Riemann) tensor 6 

absorb, via standard renormalization procedures, the divergences inherited from the point-particle 

approximation. In treating spin effects, the first higher dimensional operator we encounter is the 



self-induced term 



24 



26 



m (a = 1,2) 

L ES 2 



^ES 2 ^ab ca QC b 
2m A r^ b Ac*A 



(26) 



where E a b = R^au^v^e^v^e^ (the electric component of the Riemann tensor), and C E ^ 2 are coef- 
ficients describing the short-distance physics at the scale R. For Kerr black holes C^s 2 = 1) and 
the operator in Eq. (|26|) describes the self-induced mass quadrupole moment of a rotating black 



hole 



24, 



291 ] . In the weak gravity approximation we have at linear order in h 



fj,f 



bs 



L< 



ES 2 



L 



ES 2 



TV 2 

^ES 2 



Ces 2 
Amrrip 

Ces 2 
2mm f 

c es 2 
2mm„ 



u nik njk 



h ol;ij v l S ik S jk , 



-h mj S i0 S j0 + S ik S jk (h im v l - -h lr>ij v r v l + h Hdr v l v r + -h iU a l 



Ok ojki 



hoax + S w S^h 00tlj v 



(27) 
(28) 

(29) 



with a = -¥r the particle's acceleration, and the quadratic couplings we need are 



j h 2 _ C ES 2 ik j k 

Les2 - 8^r 



hoo,ihooj + hoo : i(hijj — 2huj) + /ioo^oo,jj — 2huhoo,ij 



+ 



(30) 



Once again we split the gravitational field into potential and radiation modes, h^ = h^ u + H^, 
to generate the Feynman rules some of which are depicted in Figs. 03&. 



Operators with Ricci tensor or scalar are redundant and can be removed via field redefinitions [23 . [52 



13 




b) 

FIG. 3: O(S^) worldline coupling to: a) radiation and b) potential modes. A box represents couplings 
quadratic in spin due to the Unite size operator in L E gi. 

D. Radiation with spin 

For the spin induced multipole moments the task is to compute all Feynman diagrams with an 
external radiated graviton and worldline spin couplings which contribute to the (spin-dependent) 
stress tensor T^(£, k) in Eq. ([3|) at the desired order. From there we obtain the multipole moments 
using Eqs. P-fTTD and (fT3|) . 

Some of these diagrams are such that no potential modes are involved, and simply consist of 
a single radiation mode coupling to each worldline independently (depicted in Figs. W& and Eh for 
the spin and finite size couplings from Eqs. (|16l - ll9p and (|27l -l29 p respectively.) We refer to these 
contributions as 'worldline radiation' since one can simply read off T^ v from the worldline action. 

The remaining diagrams represent nonlinear gravitational effects and stem from two different 
sources: 1. The explicit nonlinear terms from the spin couplings such as SHh and S 2 Hh in Eqs. 
(|20p and (|30p . 2. The bulk action, which starts with the term HHh, namely the 'three-graviton' 
vertex 7 . This is all we need in this paper. In these cases the radiation is emitted off the 
stress-energy of the system that includes the gravitational component. 

Throughout our computations we encounter contributions linear and quadratic in spin. Com- 
ponents of cubic order in spin are omitted since they first contribute at 3.5PN. We extract the 
multipole moments from moments of T^ v which in turn are computed from Feynman diagrams 
yielding T' u/ (£, k). To clarify at which order spin effects start to conribute to the multipole mo- 
ments, it is useful to know the scaling of the leading effects in T^ v . They follow from the linearized 
couplings in Eqs. (J16M18P and (|27H29|) as well as their spinless counterparts. Since even the leading 
terms in T^ u itself do not scale uniformly, we instead give the scalings for the leading contributions 
to general moments Kf" = J cPxT^x? 1 . . . x l£ in Table HI The scalings in Table U also clarify which 
moments of T^ v in the expressions for the multipole moments in Eqs. pillip first contribute at a 



7 We obtain its Feynman rule by expanding the Einstein-Hilbert action plus gauge fixing term [22 
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0(J8) 


0(S A ) 


G(S\) 


K™ 


mr i. 


mr v 


f 4 

mr v 


Kf 


mr v 


f i 
mr v 


mr v 


Kf 


i i 
mr v 


mr v 


£ 6 

mr v 



TABLE I: Scalings for leading terms in moments of T /il/ (i, k) for various orders in spin. These scalings are 
valid for all moments £ > 2. 



given order, where we remind the reader that time derivatives scale as v/r. 



Below we review in more detail which are the the necessary ingredients for the energy flux to 
3PN order. 

1. Spin- orbit radiation 



Spin-orbit effects in the energy flux at LO arise both from mass quadrupole and current 
quadrupole radiation, and give 1.5PN corrections to the LO spinless energy flux. 

More precisely, the leading 0(Sa) effects in the energy flux (at 1.5PN) can enter in three 
different ways: 1. Due to the interference of the leading spinless mass quadrupole (OPN) with the 
leading 0(Sa) mass quadrupole (1.5PN), 2. Due to the interference of the leading spinless current 
quadrupole (OPN) with the leading 0(Sa) current quadrupole (0.5PN) or 3. Due to the spinless 
mass quadrupole radiation (OPN) but using the LO spin-orbit equations of motion (1.5PN). 

Now for the 0(Sa) effects in the energy flux at NLO or 2.5PN, the additional factor of v 2 may 
enter in many ways. However, it is clear that we need the mass quadrupole and current quadrupole, 
both spinless and 0(Sa), to NLO, as well as the spinless and spin-orbit equations of motion also 



to NLO 
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30l | . Moreover, we have to include the LO mass octupole and the LO current octupole, 



both spinless and 0(Sa)- As we shall see, the matching calculation of the mass quadrupole and 
current quadrupole at NLO requires the inclusion of nonlinear effects. 

Since the leading 0(Sa) radiation effects arise at 1.5PN, and the leading tail effect yields a 
1.5PN correction to the flux, we also need to include the tail contribution to the energy flux at 



3PN linear in spin, see 



37] for a generic expression to compute such effects. 
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2. Spin-spin radiation 

The LO spin-spin contributions to the energy flux first appear at 2PN. They can for example 
arise when squaring the leading 0(Sa) current quadrupole, which gives both 0(Sa^b) and 0(S A ) 
terms in the energy flux. Note that the leading 0(Sa) mass quadrupole (1.5PN) does not contribute 
to spin-spin effects in the energy flux until 3PN. Therefore, the only 0(Sa) multipole moment 
needed to NLO for the spin-spin effects in the 3PN energy flux is the current quadrupole, whereas 
the leading 0(Sa) expressions for the mass quadrupole and the current octupole are sufficient. 

Besides the 0(Sa) contributions to the multipole moments discussed above we also have 0(S A ) 
as well as 0(S a Sb) components. The LO contribution to the 0(S A ) energy flux (at 2PN) from the 
0(S A ) multipole moments arises due to the leading 0(S A ) contribution to the mass quadrupole 
(2PN). For the 0(S A ) power at 3PN we compute the 0(S A ) mass quadrupole to NLO as well as 
the leading current quadrupole and the leading mass octupole at 0(S A ), along with their spinless 
expressions. 

In the EFT that lives between the scales R and r there cannot be any worldline coupling of 
(D(SaSb)- Thus, (D(SaSb) worldline radiation cannot exist and the Feynman diagrams for every 
(D(SaSb) component of the multipole moments must involve potential graviton lines connecting 
the two worldlines, where on each worldline there has to be a coupling linear in spin. The leading 
contribution of this sort arises for the mass quadrupole moment at 3PN. For the power at 3PN, 
we only need to calculate this leading (D(SaSb) component for the mass quadrupole. However, 
there is another 3PN contribution to the mass quadrupole at 0(Sa&b) which arises indirectly 
due to the replacement of S l ° using Eq. (|25l) in the leading mass quadrupole at 0(Sa)- This is 
similar to what happened for the NLO conservative spin(l)spin(2) dynamics where the leading 



order spin-orbit diagrams contributed a term through the replacement of S A [27 



28]. 



Furthermore, we recall that the spin-dependence in the energy flux to 3PN may always enter 



from spinless multiple moments 
derive from the potentials in 



jy u sing the spin-dependent conservative equations of motion that 
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2Sj-l3fl] to NLO. 



III. SPIN CONTRIBUTIONS TO THE SOURCE MULTIPOLE MOMENTS 

Now we set out to compute the necessary Feynman diagrams which yield all spin contributions 
to the source multipole moments for the energy flux to 3PN. The spin-independent multipoles 
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which enter the spin-dependent part of the energy flux to 3PN are calculated partly in [37f] and in 
Appendix iBl 

A. Mass quadrupole moment P' J 

We begin with the calculation of the mass quadrupole which is the most elaborate one. We 
present first the contributions from worldline radiation which arises from diagrams without any 
potential modes and then calculate the more complicated diagrams involving potential modes and 
nonlinearities. 

1. 0(S A ) worldline radiation 

The Feynman diagram which gives rise to worldline radiation linear in spin is depicted in 
Fig. |5Jz. To compute these contributions we can either consider T^ u in coordinate space or in 
mixed coordinate-momentum space. From Eq. (|15p we obtain the compact expression for the 
contribution to T^ u from Fig. [2b. in coordinate space [241 ] 

T i>) = \ E / dt ' a ^ 4 (* - x ^')) {sTit'wt') + sr (0*4(0) , (3i) 

A J 

where x = (t,x), d a = -J^ and x A {t') = (£'> x a(0)- From now on we will suppress the implicit 
time dependence of the worldline variables. The expressions in mixed coordinate-momentum space 
read 

2jg(i, k) = J2 S °A ik'e-*** (32) 

A 

T &^ k ) = \ E { S a & + S a^a M + S° A l zk • v A - S° l ) e~^ A (33) 

A 
W*> k ) = \ E { ( 5 M + S^a) *k< + jk • v A (S«vi + S / V A ) 

A 

- fiSfvi - S A W A - Sf a \ - Sjai} e- ik ^. (34) 

For the 0(S A ) mass quadrupole at NLO the subleading terms in the components of T^ v are required 
in the pieces of P 3 of Eqs. dU-[6]) which give the leading 0(S A ) mass quadrupole, 

d 3 x (t 00 + T u - -T^iA [xV] S tf, (35) 



whereas the higher order terms 
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d 3 X ( lir 00 X 2 + -f U X 2 + \f lm ^ m - V'x'x 2 ) [X'X^STF 

1 42 21 6 7 ' 



(36) 



only require the leading terms in the components of T^ v . 

Using either form of T^ v it is straightforward to compute the moments in Eqs. 
NLO we get for the mass quadrupole contributions from Fig. ^Bp, 



. Up to 



E 



2k + jH (y A x S A ) l ^ A -|(x A x S A )W j A -|(x A x Sa)V a 



■^ <{ ^t v a • xa(va x Sa) 1 x^ 



+ 



d Uk 
di\!> 

d 2 ri + iik 2 



dt 2 



v^2G 

+ K> 



.. | x A (v A x S A )V A + -x a (Sa x xaY^a 

\ . - 1 — 6« 

+ -^A ■ va(xa x SaT^a h 42~ ( Vyi X Sa ^ ' Xyl X ^ X A 

2ms ((va - vb) x Sa) 4 x^ + -=■ ( Sa • rS#x A - Sa • S B r l x 



STF 



AS 



STF 



+ 2k(1-k) J^ 



^(vaxSa)V, 



(37) 



STF 



where we applied the replacement rule for S A in Eq. (|25p and used S^ = e* J S A . The corrections 
due to Eq. (|25l) are the terms in the last two lines of Eq. (|37|) , and we notice that this replacement 
also induces an O(SaSb) contribution to the mass quadrupole which is of order 3PN and which 
we have to keep. Furthermore, note that we kept a term proportional to Sa- Upon use of the 
spin equations of motion (which we will not perform until we compute the energy flux in 31]) this 



term becomes an 0(Sa) contribution at 2.5PN and spin-spin contributions at 3PN to the mass 
quadrupole moment. Employing the covariant SSC k = 1, Eq. (|37|) becomes 



% ~ E 



8 ■ • 4 4 4 d 

-(va x Sa)'x^ - -(xa x S A ) l v A - -(xa x Sa)V 



'^ " \jt l VA ' XA(VA X SA)i ^} 



7dt 2 



1 



+ ^~n?. <i 4x a( v a x Sa)*x a + ^xa • v A ( Xj4 x Sa)Sc^ 



+ x a (Sa x x A )V A - -(va x Sa) • xa x a x 



STF 






A,B 



2m B ((va - v B ) x Sa) 4 x a + -=■ ( Sa • rS#x A - Sa • Ssr r x 



(38) 



STF 
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2. 0(S\) worldline radiation 

The Feynman diagram for 0(S\) worldline radiation is shown in Fig. [3b- The relevant couplings 
arise due to Eqs. ([27I - I29p . The expressions for T^ v from Fig. [3b are 



2f(t,k) 



E 



c 



£S 2 



2rriA 



1 + Li ) £^ fc _ 5 ;5 / - 25'MM 



M ^A ^A J A 



A°A V A 



k i k i e -*.x^ 



c 



(^4) 



E ES 2 qi l qkl i \.j\M -iV--x. A 



2m A 



E 



c 



ES 2 



2l7lA 



okl qkl -l, „ , nWomL^iium , nqkmqlmA -\,k 
-&A >^>A l * ' a A + ^A '-'A v A^ **■ + ^A '-'A a A z * c 



-ik-X4 



(39) 
(40) 
(41) 



where we dropped time derivatives of the spin variables which do not contribute at the order we 
are working. From here we obtain the moments 



^f [ X ^] S TF=E 



C 



(A) 
ES 2 



2mA 



_ 2 + {2k 2 - l)v A ) S A S A - 4k(k - 1)S A • v A S A v A 
+2k(k - 2)S A v A v A 



STF 



fyjfa [xV] 

d 3 xr^x' [x^] 



STF 



STF 



E 



C 



ES 2 



m A 



-viS^-2S A -a A S A 4 



A 

" y^ °es 2 
<i£ t-* 1 mA 

A 



STF 



(42) 
(43) 



x A • v A Sf Sj + 2S A 'v A ^x A " 



STF 



E 



c 



(A) 

ES 2 



m A 



-(vi + x A • a A )S A S A + 2S A (v A v A + x A a A ) 
-2S A • v A S A v A - 2S A • a A S A x A " 



^«- 2 [^]stf = 4E 



C 



{A) r 
ES 2 



STF dt 2 '-? mA 

A 



S4^A ^A 



STF 



kl k qil^j _i_ okl okl i j 



(44) 



2 OtkqJK , a qkl k qll 3 _|_ qkl qkl I 



A X A^A^A ~ ^A ^A -V1 A A 



STF 



E 



C 



(A) 
ES 2 



m A 



-2(v A + x A • a A )S A S J A + 12Si(vXt + *W) 



-8S A • v A S A v A - 4S A • a A S A x J A - 4S A • x A S A a A 



STF 



(45) 
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which we combine to the mass quadrupole using Eqs. (jU [5]). We obtain 



& 



E 



c 



(A) 
ES 2 



ruA 



v J .\ \ - J + V A ( « - 42 J + 2i aA ' XA J ~ T7f X A S A a A • Sa 



22 



29 



17 



21 



+ S 2 , 



and in the covariant SSC we have 



^a A S A S A • x A + I - - 2(k - 1)k ) v A S A S A • v A 
10 \ , ,• 10 



21 



k(k - 2) + ^ ] v A v A + — a A x A 



21 



STF 



4 - s 



c 



(A) 

ES 2 



m A 



S A S A 



13 2 17 

42 VA+ 21 ; 



I + — v i + §j a A • x ^ 1 + S A ( - — v A v J A + — a A x J 
4 , _,- „ 22 



1 I .• j 10 
21 A A 21 



^-x^S^aA • Sa + ^v A S A S A • v A - ^a A S A SA ■ x A 



STF 



(46) 



(47) 



3. Nonlinear gravitational contributions at 0{Sa^>b) 

Let us start the discussion of nonlinear diagrams with the new results quadratic in spin and 
first consider the spin(l)spin(2) contributions. For the leading 0(SaSb) mass quadrupole moment 
(which is of order 3PN) we need to calculate Feynman diagrams which involve the exchange of a 
potential graviton shown in Figs. 0Jz6. 





S'2 

a) 



FIG. 4: Nonlinear gravitational effects at 0(SiS2)- 



The radiation graviton can then either couple to the worldline together with a potential 
graviton (Fig. 0b.) or to the potential graviton away from the worldline (Fig. 0b). 



Let us calculate in somewhat more detail the contribution to I lJ from the diagram in Fig. \$o,, 
beginning with T 00 . For this we need the HSIiqq coupling in Eq. (f20|) . concretely 4 -^2-S , f-ffoi^oo,Z) 
that we ought to contract with the LO spin vertex in Eq. ()16j) . Then the contribution to the 
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effective action in Eq. ([2D will be of the form (in momentum space) 

- ^ J d 3 kT^(t,k)h 00 (t,-k) = _i f -!-5fS 2 fcm (F 0i (x 1 )Fo fc , m (x 2 ))/ i oo,Kxi). (48) 
Using 

<flbi(x A )2r fc,j(*B)> = P 0t ok I ^3 ^ fa') e-*i-^-^) = i^^™, (49) 

we have (after summing over mirror images) 

T i( k ) = E S « r x S «) x S ^ Hk')e- k ^. (50) 

A,B r 



Similarly, for T we need the coupling -t—^SHq^ j/i/& from Eq. (|20p . Following the same steps 



p 



the result is (after tracing over Ik) 



r A( k ) = — r E ( s ^ • s ^ - 3S ^ • nS ^ • n ) e " 4kXA - ( 5l ) 



r A,B 



Since the second and third terms in Eq. ([5]) are clearly subleading, to compute the contribution 
to the mass quadrupole we just need to expand Tr? (k) + 2m(k) to second order in k and read off 
the contribution to the mass quadrupole using (fT3|) . The result is 



2(7 r . • 1 

J^, = y^— Sa ■ rS B x A - S B • Sa^x^ - (Sa • S B - 3Sa ■ nS B • n) x A x A . (52) 
A,_B 

For the diagram in Fig. ®> we need the three-graviton coupling as well as two LO spin insertions 
of Eq. (I16p . This diagram is more symmetric, and since we only need T 00 + T M , we can combine 
both Feynman rules, i.e. we add the one with an external /loo an d the one with an external ha- 
which we trace over kl. The expression for the combined three-graviton vertex contracted with the 
two potential propagators reads 

(H 0i (p + k)H oj (p)[HHh(k)}) = ^-[2kV + 2(pV - kV)]-^— =^, (53) 

n 

where [HHh(k)] represents the three-graviton vertex [221 ] . The open indices are contracted with 
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the spin couplings on each worldline, then we have 
(T 00 +T k %(k) - 



rl 3 r> 1 1 1 

1 vkH(p ! y-k 4 P j )]- 



(2tt) 3 4m 2 , 



V(k + p) 2 



,SiV^'(-P Z -k')e ip - Xl e^ (p+k) - X2 . 

(54) 



Expanding to to second order in momentum and after some cancelations we have 



(T uu + T ^ (k) 



d 3 P ' gikgjl 



(2vr) 3 2ml 
and performing the Fourier integration we obtain 



k'k J p p ■ Xi ipX2 
P 4 



+ ..., 



AG 



l I ~,2 XVni „m_J 



(T uu + r* fc ]^(k) = -jSrS 3 2 L (r 2 5 ml - r m r 



yy 



H • 



Therefore from Eq. (fT3 






STF 



A,B 



The combined result from these diagrams reads 



G 






12x^S^ • nSfl • n - 4S,! • Sb^b + 4 Sa ■ r(2S' B xj, - x^S^ 



(55) 



(56) 



(57) 



A,B 



STF 



(58) 



^. Nonlinear gravitational contributions at 0(S A ) 

We now calculate the finite size contributions. Even though the leading 0(S A ) mass quadrupole 
is entirely due to T 00 worldline radiation, we need to consider nonlinear contributions not only to 
T 00 but also to T %3 . This is because at O(S^), the leading T l i is not purely due to worldline 
radiation but gets a contribution from diagrams with potential exchange as well. This is the same 
as what happens in the spinless case 37|. The three Feynman diagrams with potential exchange 
which enter the 0(S A ) mass quadrupole at NLO are shown in Fig. [5b&c. 

The calculations of the Feynman diagrams proceed analogously to the nonlinear ones at 
0(SaSb) above. Notice that the first diagram in Fig. [5ji only contributes to T 00 since the worldline 
coupling to hij introduces an extra power of v 2 making it a higher order term in T l K We find 



^xTjgfxVjxF 



^-^ tuaC ES 2 G 



A,B 



m B 



2r 3 



(S| - 3(S B • n) 2 ) [x^] TF 



(59) 




a) 



C, 



Ci 



rii-z 




111-2 




FIG. 5: Nonlinear gravitational effects from finite size at 0(S A ). 
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d *% [x x^]tf - 2^ 9 _3 



A,B 



Itxibt 



4S|x^ B + 2r 2 S^S J B - 3 (S| + (S B • n) 2 ) x*,4 (60) 



STF 



d 3 x2[|VxV 



£5 2rn ^ 



ft 2m B r3 



4S B • rS^ - 2 (S| - 3(S B • n) 2 ) x^ 



STF 



(61) 



and 



^3 x(T oo + T , fcfe[x , xV = ^ ^ [-2(S|-3(S B -n) 2 )(x^ + x A x A ) 

3 A,B Br 



+4r z S B S J B - 8Si,r l x^ + 8S B • rS^ 



STF 



(62) 



For the total contribution to the 0(S A ) mass quadrupole from Figs. \5jp,bc we have 



K7 = V^ GC ES2 m A r 

"ISJlfec Z— • 9-m n r3 



A,B 



2mBr- 



(S| + 9(S S • n) 2 )x^ + (3(S B • n) 2 - S|) x^ 



+6r z S J B S J B - 4S^x A x^ + 12S B • rS^ 



STF 



(63) 



5. Nonlinear gravitational contributions at 0{Sa) 



Finally for the nonlinear 0(Sa) contributions there exist three Feynman diagrams as seen in 
Fig. [6)260. 

For the NLO mass quadrupole at 0(Sa) we also have to calculate subleading terms for the 
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FIG. 6: Nonlinear gravitational effects at 0(Sa)- 



moments in Eq. f|35f) . i.e. for all components of T^ v '. This means that we need to consider the 
diagrams in Fig. [6] for external radiation field legs with /too , ^Oi an d h%j ■ The diagram in Fig. 
02 only contributes at the order we are working for Tiqq since the other components introduce an 
additional power of v 2 in the couplings to the worldlines, hence 



X _B X B 



() , d 3 x 3^[xV]tf = - J2 —IT- K 1 + K )( v ^ >< S a) • r - 2(v B x S A ) • r} 

J A,B T 

For the contributions from the diagram in Fig. (6f> we have 

f d 3 x (r 00 + r-)^ [xV] TF = J2 = ^r- [i 2 ^ x s ^) ' r - A ^ B x s ^) • r > x a4 

J A,B r 



(64) 



(« 



-4r 2 (v A x S A ) 4 x^ + (v fl x S^M 



STF 



(65) 



and 



^^TJSv'fJJi ° V^ GrriB 



dt ^ r 3 
AB 



-2r 2 (xA x SaT^a --(xjjx Sa) • xax a x^ 



• (66) 



J STF 
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Finally Fig. [Ob yields 



A,B 



4G7TT-B 



(67) 



(v B x S A )*(x A + x B ) + (1 - k)(v a x S A )V A + (r x S A y (V B - v A - ^(4 + 4) 



1 



+ 2^( VA X Sa) " r { (1 + K ) x «4 + ( K - !) x Axi} 



STF 



and 



d ^— v Gm^ 



7 " E A,B 



r 2 -(r x S A )V + (x fl x S A )V A + x J B ) (68) 



(x A x S A ) • x B x B x^ - r • x B (r x S A ) J (x A + x B ) 



STF 



Combining these ingredients we have in the covariant SSC 






2Gm B r 



A,B 



(v B x S A ) • r(x B x J B - 2x A x A ) + (v A x S A ) • r(x A x A + x B x J B ) 



+ 2;- \ (vb x S A )% - (v A X S A )'x A + (r x 8 A )* (v£ - v A - ^(x A + x s )) } 



STF 



4y^ d 

3 , ^ 
A,B 



Gms 



, 1 r 2 ( ^(xb x S A )'x A -^(x a x S a )*x a + ^(xb x S A )S4, -^(x A x Sa)*xJ, 



-r • xb (r x Sa)^x a + x^) + -(x A x S A ) ■ x B (x A x A - 2x i B x' B ] 



(69) 



STF 



B. Current Quadrupole 



1. Worldline radiation 



The current quadrupole J* J at 0(S A ) is needed to NLO. The components of T^ v at 0(S A ) from 
worldline radiation including all subleading terms were already given above in Eqs. (|3ip or (J32I - 
34"|) . Expanding these to obtain the moments of T^ v in the current quadrupole J y of Eqs. ([9lll0p. 
we find for the worldline contributions to the current quadrupole from the diagram in Fig. [2k 



% - J2 



14 



1 



S A • x A a A x A + -x A x A a A • S A + I — - - J S A • v A v A x A 



1 K 



+ ^S>M + Si4Q-(^-f)v| 



K 3 

2+ n )^-x A 



+ - - « 1 v A S A v A • x A + I ac - — S A • x A v A v A 



14 



(70) 



STF 
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which takes the following form in the covariant SSC 



% = Yl s ^ x a U + fa - 7 aA ' x ^) " 7 v ^ s >a • xa + §s» A 



11 



+-Sa • x A a A x A + -x A x A aA ■ S A ~ =Sa ■ vav a x a + — Sa • x A v A v 



'^ 



(71) 



STF 



The worldline radiation contribution from the diagram in Fig. [3^ to the current quadrupole J lJ 
at 0(S A ) is just needed at LO and follows from Eqs. (|4"0]) . This yields 



%-z2 



C 



(A) 
ES 2 



m A 



( VA x s A ys 3 A 



STF 



(72) 



£. Nonlinear gravitational effects 

For the nonlinear gravitational contributions to the current quadrupole we only need to compute 
the diagrams in Figs. 0)c, since the contribution from [6b can be shown to be subleading. In order 
to calculate the T° l that enters in Eq. ([9|) we are only interested in diagrams with an external JiQi. 
Notice that we already performed this computation for the NLO contributions linear in spin for 
the mass quadrupole moment. Thus the results for the current quadrupole from these diagrams 
are 



and 



£ 

A,B 



GrtiB 



S A • x A rV A - x A S J A (6r 2 + x A • r) 



STF 



(73) 



% - J2 



A,B 



GrriB 
2r 3 



-3Sa • x B (x A x A - x B x B ) + S A • x A (x A x A + 2x A x B - 3x B x B ) 



+2S A x A x A • r 



STF 



(74) 



which add up to 



s« 



EGrriB 
A,B 



Sa • x a ( 2x A x A + x A x B — 3x B x B 



+3S A • x B (x B x B - x A x A ) + S A x A (x A • r - 6r 2 ) 



STF 



(75) 
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C. Mass octupole 

For the mass octupole we only need the pieces linear in h^ v coming from the worldline couplings. 
Using Eqs. (0 [SJ) and (J3TJ) we obtain 



« = E 



I + 3k ) (v A x Sa) V A x£ - 3(x A x S A )\ J A ^ A 



(76) 



STF 



and 



f 



3C 



(*) 



^^Mi [S ^S^x fc . 



m A 



a°a x aJstf- 



(77) 



D. Current octupole 



Similarly to the mass octupole, here we only need the linearized worldline coupling. Using Eqs. 
TTT) and (I3TT) we obtain 



-%, ~ 2 Z> [ X A X A S A 



STF 



(78) 
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IV. TOWARDS SPIN EFFECTS IN THE ENERGY FLUX TO 3PN 



Let us start by collecting the different contributions to the multipole moments in Eqs. (|38l 1471 
EH Efl E3 EH EHD. In the covariant SSC the results are 



J l 3 






-(v A x Sa) j x^ - -(x a x S A )V A --(x a x S a )*x 



« V J 



(79) 



r ■ ■ -I 1 d 2 f 1 

|v A • x A (v A x S A ) J x A j + ---2 < -x A • v A (x A x S A ) l x A 



+ 4x 2 4 (v A x SaYx-a + x A (S A x x A )V A - ^(v A x S A ) ■ x A x A x^ 



+£ 



2Gm B 



STF 

(vfl x S A ) • r(x J B x J B - 2x A x^) + (v A x S A ) ■ r(x A x^ + x^x^) 



A,B 



+2r 2 {(v B x S a )Vb - x J A ) + (r x S A )* (vj, - vj - ^(4 + 4)) } 



STF 



2^ d 

3f^dt 



Grris 



{r 2 ((x B x S A )V A - 3(x A x S A )V A + 3(x B x S A )% - (x A x S A ) V B 



-2r • x B (r x S A )*(x A + x J B ) + (x A x S A ) • x B (x A x A - 2x 4 B x^)} 



+£ 



C 



(A) 

ES 2 



m A 



&A&A ( - 1 + ^A + ^a A • X A 1 + S A [-^>A + ^ a W A 



32 f ll , T i < 



STF 

10 



q 4 • 22 

— x A S^a A • S A + -v A S^S A • v A - — a A S^S A • x A 



+ Z-> O.r-3 



AS 



^^ (8% + 9(S B • n) 2 ) x^4 + 6 ^^^-'« «■' 



rrt-B 



m B 



STF 



+ 



ES2 A (3(S B • n) 2 - S|) + 12S A ■ nS B • n - 4S A • S B x A x J A 
m B J 

C ES2 m A 

TUB ..-.,, 

STF 



S|x A x J B + 4 |3 ^g^^ S B • r + 2S A • r | S B ^ B 
\ mB / 
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E 






( Vj4 x Sa) j S A + S A x A I - + -v A - -a A • xa j - -v A S A v A • x A 



11 2 ■ ■ 1 • ■ 3 -11 

+^ s >a x a + ^ S A • x A a A x^ + -x A x^a A • S A - -S A ■ vav a x a + — S A • xav a v a 

^-\ Cms 



2r 3 



3S A • x jB (x i B x^ - x A x A ) + Sa • xa(2x a x a + x A x B - 3x s x s 



+S A x A (xA-r-6r 2 ; 



J, 



i)k 



S A 



rfjfc 

'S A ,S2 



; EK* 



STF 






STF 



E 



^) 



^(va x SA) i x A x A - 3(x A x Sa)V a x a - 3^S A S A x A 



STF 



(80) 



STF 



(81) 
(82) 



plus the spin-independent contributions (see Appendix [B]) 



vs 



4 



E m ^ 

A 



i+|vi-E^)<4+ 



Y^ m A ( 1 + 



11 d 2 r 2 i , i 4 d r 
^^|xaXaXa)-3^{xa-vax,, 



x J a} 



¥ [ 



(xa x va) 1 x a 



STF 



+E 



GrriArriB 



2(xa x v A )*x A 



(83) 
TF 

(84) 



A,B 



11 



— (x B X V A ) l X B - -(xb X V A ) 4 X A + (x A X V A )Vg + -(x A X X B )V 



»„J 



4 
va r 



4r^ 



(xa x xb) 1 (x a + x b ) 



1 d 

stf + 28 dt 



^ wa(xa x va)*(3x a v a - xa • vax a ) 



GmArriB 



+ E — Sla — x a (xa x xb) j (6x a - 7xa ■ x B + 7x. B) 



A,B 



2r 3 



STF 



J, 



«J& 






^2 m A [(XA X V A ) l X A X A 
A 

Y, mA 

A 



STF 



X A X A X A 



TF 



(85) 
(86) 



The spin dependence of the multipoles is given in the covariant SSC which was used in [28H3C | 
to compute the NLO spin(l)spin(2), spin(l)spin(l) and spin-orbit potentials. We could combine 
the results also in the NW SSC, k = 1/2. This would simplify the calculation of the power loss 



2fi^y,yi. 



However, the NW SSC 



for spin(l)spin(2) and spin-orbit effects by using the results in 

ceases to be useful once we include spin(l)spin(l) (finite size) corrections, where we found the 

covariant SSC to be more convenient. In any case, we can always transform the spin dynamics into 



precession form at the end of the day by means of a spin (and coordinate) redefinition 



22 



28M30I]. 



Ultimately the power loss will not depend on the choice of coordinates but its specific form will 



29 



depend on the choice of spin variables. 

Our results in Eqs. (|79l - l86j) are the last missing ingredients for the spin-dependent parts of 
the energy flux to 3PN. Its computation will be presented in |3ll |. Besides the multipole moments 



in Eqs. (|79H86p we also need to make use of a) the conservative spin dynamics to NLO 24|, I26H29 ] 
where appropriate, and b) the modification of the energy flux linear in spin due to the leading tail 
effect [37]. Furthermore, we can also use the multipole moments presented here in the computation 
of the gravitational waveform hj- . 



V. CONCLUSIONS 



In this paper we computed the necessary multipole moments, expressed 



the covariant SSC, which together with the results reported in 



22 



24 



in Eqs. (|79l - ISo]) in 



28 30, 



37J will lead to the 



spin-orbit and spin-spin contributions to the power loss to 3PN order via standard procedures. For 



,y,r 



the computation we employed the EFT techniques developed in 22j, [24], [37]]. There is clearly no 
stumbling block to continue this program to higher orders in the PN expansion. We will report 
the final expressions for the energy flux and the phase in a separate paper [311 ] . 
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Appendix A: Contribution from the term responsible for the preservation of the SSC 

describe spin dynamics in GR, and ensure the preservation of the SSC upon evolution, in 



25, 



28 



291 ] the following Routhian was introduced (up to finite size effects) 



n = ~ E ( ™>A\fA + IsfuJa^u'X + J-R deab ( XA )S Ac d Sf^^ + •••), (Al) 

where the ellipses represent nonlinear terms in the curvature necessary to account for the mismatch 
between p and u once the SSC is enforced. Since at 3PN order we can consider the covariant SSC 



30 



to be S ab Ub = 0, we can ignore the higher order terms. The equations of motion follow from 

6 (' .. . dS ab 



Sx -j KdX = ' -*r = ( s •*>• (A2) 

and the Riemann dependent term in the Routhian guarantees the SSC is preserved upon evolution 



28, 



29]. Notice that it vanishes once the SSC is enforced, however, it gives a nonzero contribution 
to the equations of motion from Eqs. ()A2p , In other words, we cannot set the SSC to zero until 
after we perform the variation of the Routhian. In particular, 

P a = -L= ( mu a + -LRp upa S a PSP°uA . (A3) 

Naively it appears as if this extra term in the Routhian does not contribute to the radiation 
sector, since the SSC is conserved in time. However, we know it has to contribute to the stress- 
energy tensor, from which the equations of motion derive. Once again the subtle point lies in the 
fact that we cannot apply the SSC before we expand the metric field in terms of the radiation and 
potential modes. We notice that only the term where the radiation field stems from the vierbein 
inside the SSC contributes to the one-point function (see Eq. [3|), and it reads 

^ L -R d£ab (H llv (x A ))St^=Sfh fM: u l X, (A4) 



^2m/™^^' ~*vsr 



where we used 



<% = % + & + '" ( A5 ) 



1 

— / 
2 

All terms besides the one in Eq. (|A4p which result from expanding the third term in Eq. (jAip to 
linear order in the radiation field are proportional to the SSC and do not contribute. 

The expression in (|A4|) produces a ddHh type of coupling which leads to a diagram equivalent 
to Fig. 0>. By standard power counting it is easy to show that this extra piece contributes an 
0{S A ) term in the mass quadrupole moment at 4PN. 

Notice that the coupling in (|A4j) leads to the following expression for the stress tensor (after 
adding the spinless part from the worldline) 

T ^RSSuu = Ys mAU A U A + 2^( R P*P°)( XA ') U A S A fi S A U A = J2 P A U a\ ( A6 ) 



31 



where the expectation value of the Riemann tensor, i.e. {Rpapa), is obtained in the background 
field of the companion. This is in accordance with similar expressions found in the literature (see 



for instance 



20, 



2i|- 



Appendix B: Spin-independent multipole moments 



In this Appendix we summarize (within our conventions) the spin-independent multipole 
moments necessary to compute the radiated power linear in spin to 3PN order, some of which 
have been already reported in [371 ] . 




rri\ 




m-> 



a) 



b) 



FIG. 7: Diagrams representing the spin-independent nonlinear gravitational contributions to the multipole 
moments. 



lid 2 



For the mass quadrupole (see Eq. (101) in [3J 

4ii = E m ^ f (i + |vi) x^ + -^ {xix^4} 






TF 



(Bl) 



to 1PN order. For the current quadrupole the LO term reads (see Eq. (89) in 



37|) 



■#' = £ 



m^ 



(x A x v A )*x^ 



STF 



(B2) 



whereas for the NLO contributions the new ingredients follow from: the radiation off the worldline 
couplings 



/ d 3 x2j#x J x' = y~]mA— v^x^x^, 



(B3) 
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and the diagrams in Figs. \T)flb with an external h t 



Ok 



I 



txlftM = ^°^vb& A , 



(B4) 



A,B 



d 3 xT^x j x l 



2Gm A rriB 



E 

A,B 



E 

A,B 

GrriArriB ' 



v k A x 3 B x B + v k w\(x A + x B ) 1 



(B5) 



Qr 



Va(^x^ + x J A x A + x> A x l B ) + r k [2v j A x l A + v^x 



^A ' *" / j i j I j I ' 

(x A x A + x B x B + x A x B 



We also need the contribution to the current quadrupole from Eq. (|10p with the LO T y . The 



LO T lJ is comprised of a worldline contribution 



m 



TV(t,k) = y / m AVA v j A e- ik ^, 



A 



(B6) 



as well as of a term from the diagram in Fig. [7f> with an external hij . It yields the moments of T lJ 
for Eq. (PBJ 



1 

28 



E «fc j d z xf^x j x k x r ' 



STF 



28^^^ £ [ (XA X Xfl) '^ ( " XA " XB + X|) 



AS 



STF 



(B7) 



and 

3_ 

28 



; «fc /, 



e ilk I d 3 xlgx l x 2 



STF 



1 d Gm A mB v-^ r. 



x A x x B )V, (6x A + 6x| - 6xa • x s ) 



STF 

(B8) 



so that the total contribution from Eq. (|l(jp reads 



J? 



dt 



^2^-( x ax v A ) l (3x A v A - x A • v A x A ) 



28 



+ E !!! S^^(x A xx^(6xi-7x A .x B + 7x|) 



A,B 



56r 3 



(B9) 



STF 
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Adding the pieces together to 1PN order 



4i-i = E m A ( i + ^ ) [ (xa x VA M] + £ 



GraATfiB 



STF 



A,B 



2(x A x Vi ) ! x^ (BIO) 



11 3 ■ • 7 

- — (x B x v A ) ! x J B - -(x B x vaT^-a + (xa x va) 4 x^ + -(x A x x B )V^ 



+^( x ^ xx £)Va+ x b: 



1 d 

STF + 28 dt 



^]myl(x4 x va) 1 (3x a v^ - xa • vax-^) 



+ E ^StM( x * x xb)^(6x a - 7xa • xb + 7x|) 



A,B 



2r 3 



STF 



Notice that the last term in Eq. ()B10p vanishes in the CM frame and therefore does not radiate. 
Finally for the octupole moments at LO 



T ijk 



E 

A 



m A 



x A x J A x A 



TF 



(Bll) 



and 



J o k = E mA \( XA x vaT^A 



STF 



(B12) 



As a check, one can show that Eq. (IBlOj) reproduces the result in 8J for the radiated power 
that follows from the current quadrupole piece in Eq. ( [2"]) . 
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